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ORBITAL DEPENDENCE OF GALAXY PROPERTIES IN SATELLITE SYSTEMS OF GALAXIES 
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ABSTRACT 

We study the dependence of satellite galaxy properties on the distance to the host galaxy and the 
orbital motion (prograde and retrograde orbits) using the Sloan Digital Sky Survey (SDSS) data. 
From SDSS Data Release 7 we find 3515 isolated satellite systems of galaxies at z < 0.03 that contain 
8904 satellite galaxies. Using this sample we construct a catalog of 635 satellites associated with 
215 host galaxies whose spin directions are determined by our inspection of the SDSS color images 
and/or by spectroscopic observations in the literature. We divide satellite galaxies into prograde and 
retrograde orbit subsamples depending on their orbital motion respect to the spin direction of the 
host. We find that the number of galaxies in prograde orbit is nearly equal to that of retrograde orbit 
galaxies: the fraction of satellites in prograde orbit is 50 ± 2%. The velocity distribution of satellites 
with respect to their hosts is found almost symmetric: the median bulk rotation of satellites is — 1 ± 8 
km s~^ . It is found that the radial distribution of early- type satellites in prograde orbit is strongly 
concentrated toward the host while that of retrograde ones shows much less concentration. We also 
find the orbital speed of late-type satellites in prograde orbit increases as the projected distance to the 
host (i?) decreases while the speed decreases for those in retrograde orbit. At R less than 0.1 times 
the host virial radius {R < 0.1rvir,host) the orbital speed decreases in both prograde and retrograde 
orbit cases. Prograde satellites are on average fainter than retrograde satellites for both early and late 
morphological types. The u — r color becomes redder as R decreases for both prograde and retrograde 
orbit late-type satellites. The differences between prograde and retrograde orbit satellite galaxies may 
be attributed to their different origin or the different strength of physical processes that they have 
experienced through hydrodynamic interactions with their host galaxies. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: general - galaxies: morphology 
- galaxies: properties 
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1. INTRODUCTION 

In the hierarchical picture of galaxy formation, mas- 
sive galaxies are formed through merger and accretion. 
Satellite galaxies associated with massive host galaxies, 
memorize the history of their accretion and interaction 
with host galaxies. In addition, they are distributed out 
to a few virial radii of the host galaxy. Therefore, they 
are excellent probes of the dark matter halos surrounding 
the host galaxies. 

A great deal of effort has been made to understand the 
spatial distribution of satellite galaxies around their host 
galaxies. The radial distribution of the observed satel- 
lites was found to be more concentrated than that of sub- 
halos aroun d galaxy-sized halos in N-body simulations 
(jChen et al. 2006), but consistent wi t h dark matter dis - 
tribution (Ivan den Bosch et al.l 120051: IChen et al.ll2006D . 
The studies of color dependence on the radial distribu- 
tion of satellites revealed that the radial distribution of 
red satellites is more centrally concentrated than that of 
blue satellites dSales et al.l 120071: lChenll2008D . Similarly, 
the radial distribution of early morphological type satel- 
lites is mo re central ly concen trated than that of late-type 
satellites (|Ann et a L 2008; Wang et al.ll2010f ). The angu- 
lar distribution of satellite galaxies is a.lso an issue, which 
has been extensively addressed (e.g., iPrada et al.l 120031 : 



Sales & Lambas 2004. 20051: [Ann et al.l2008HBailin et~all 
2008: .Wang et al...2009. ). 
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Another interesting issue is the kinematic properties of 
satellite galaxies. Satellite galaxies can be divided into 
two classes based on their orbital motion: those in pro- 
grade orbit that rotate around the host in the same direc- 
tion as the rotation of the host, and those in retrograde 
orbit that rotate in the opposite direction. Several im- 
portant questions have been raised for these two classes: 
(1) what the number ratio of galaxies in prograde orbit 
to those in retrograde orbit is; (2) what the amplitude of 
mean bulk rotation of satellites is; and (3) whether there 
is a difference in physi cal properties between them. 

IZaritskv et al.l (J1993D searched for the satellites around 
late-type spiral galaxies with a luminosity similar to the 
Milky Way, and found 69 satellites around 45 host galax- 
ies (~1.5 satellites per host). For 35 satellites whose or- 
bital velocity directions were determined, they detected 
a small amount of rotation for the satellite systems, 
~ 29 ± 21 km s~^, in the same sense as the rotation of 
hosts' disk. They also found that the number of satellites 
in prograde orbit is nearly e qual to that in r etrog rade 
orbit. However, in NGC 5084. ICarignan et al.l ()1997n de- 
tected nine satellite galaxies, and found a clear excess 
of satellites in retrograde orbit (eight in retrograde orbit 
and one in prograde orbit). Using the ex panded catalog 
of 115 satellites around 69 host galaxies, IZaritskv et al.l 
(|1997f ) found that the number of galaxies in prograde or- 
bit is almost the same as that of galaxies in retrograde 
orbit. 
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Recently, lAzzaro et al.l ()2006D measured the direction 
of rotation fo r some host galaxi es in the list of satel- 
lite systems of lPrada et alj ()2003l ) . and secured 76 satel- 
lites associated with 43 host galaxies of which orbital 
directions were determined. With this sample, they 
found that the fraction of galaxies in prograde orbit 
(/pro) is ~ 0.6, which is sl ightly larger than th e previ- 
ous results (/pro ~ 0.52 in iZaritskv et al.lll997() , but is 
broadly consistent with those from their TV-body sim- 
ulation (/pro ~ 0.55 — 0.60). In the catalog of satel- 
lite syst ems extracted fro m the Sloan Digita l Sky Survey 
CSDSS : lYork et al.ll2000l) Data Release 6 (jBailin et al.l 
I2008D . iHerbert-Fort et al.l (|2008D determined the direc- 
tion of rotation of some host galaxies by spectroscopic 
observations. They obtained the data for 78 satellites 
associated with 63 hosts. Combining the s e dat a with 
those for the satellites in 'Za ritskv et al.l ()1997l ). they 
found that the fraction of galaxies in prograde orbit is 
still 0.53 ± 0.04, and mean bulk rotation of satellites is 
37 ±3 km s~^, in the same sense as the rotation of hosts' 
disk. Interestingly, they found that the peculiar velocity 
distribution of satellites is not a single Gaussian, but is 
described as the sum of two Gaussians. 

As summarized above, most studies focused on the first 
two issues, the fraction of galaxies in prograde orbit and 
the amplitude of mean bulk rotation of satellites. The 
difference in physical properties between them, which is 
very important to understand the formation of satellite 
systems, has not been explored much. Moreover, the 
direction of rotation of host galaxy that is essential to 
know the orbit of satellites, has been determined by the 
expensive spectroscopic observations. 

In this paper, we determine the direction of rotation 
of host galaxies using the color images provided by the 
SDSS without conducting spectroscopic observations (see 
§2.31 for details) , and present the results of a study on 
the dependence of various physical properties of satellite 
galaxies on their orbits (prograde and retrograde orbits). 
In this paper we mean by 'satellites' as the galaxies that 
are much fainter than their host galaxies that are isolated 
from other bright galaxies. Therefore, they are like the 
conventional satellites associated with the Milky Way or 
Andromeda, and are rather different from the 'satellites' 
used in the N-body simulations or group/cluster studies 
where all galaxies (or dark halos) other than the bright- 
est one (or the most massive one) are called satellites. 
Section [2] describes the observational data used in this 
study. Orbital dependence of the physical parameters of 
satellites are given in [J3] Discussion and summary are 
given in 21 ^nd iJ51 respectively. 

2. OBSERVATIONAL DATA SET 
2.1. Physical Parameters of Galaxies 

We use a spectroscopic sam ple of galaxies in the SDSS 
DR7 (jAbazaiian et al.ll2009[ ). The physical parameters 
of galaxies that we consider in this study are r-band ab- 
solute Petrosian magnitude (Mr), morphology, {u — r) 
color, equivalent width of Ha emission line [W(IIa)], 
color gradient in (g — i) color, concentration index (cin), 
internal velocity dispersion (cr), and Petrosian radius in 
i-band (i?pet)- Here we give a brief description of these 
parameters. 

The r-band absolute magnitude Mr was computed us- 



ing the formula, 

Mr = m, - 51ogio[r-(z)(l + z)] - K{z) + E{z), (1) 

where r{z) is a comoving distance at redshift z, K{z) 
is X-correction, and E{z) is the luminosity evolution 
correction. We adopt a flat ACDM cosmology with 
density parameters J7a = 0.73 and Vim — 0.27. The 
rest-frame absolute magnitudes of individual galaxies are 
computed in fixed bandpasses, shifted i o z = . 1, us- 
ing Galactic reddening correction (Schl eeel et al.lll998fl 
and X-corrections as described by Blanton et al.' ('2003') . 
The evolution correction given by Tcgmark ct al. (2004), 
E{z) = 1.6(2; — 0.1), is also applied. The distance r{z) 
has a unit of /i~^Mpc and the corresponding b\ogh in Mr 
will be omitted in this paper. 

We adopt the galaxy morphology information from the 
SDSS DR7 release of the Korea Institute for Advanced 
Study value-added galaxy catalog (KIAS-VAGC; Choi 
2010, in prep.). In this catalog galaxies are divided 
into early (ellipticals and lenticulars) and late (spirals 
and irregulars) morphological types based on their loca- 
tions in the {u — r) color versus {g — i) color gradient 
space and also in the i-band concentration index space 
(jPark fc Ch oi 2005). The resulting morphological clas- 
sification has completeness and reliability reaching 90%. 
We perform an additional visual check of the color images 
of galaxies to correct misclassifications by the automated 
scheme. In this procedure we changed the types of the 
blended or merging galaxies, blue but elliptical-shaped 
galaxies, and dusty edge-on spirals. 

The '^'^{u — r) color was computed using the extinc- 
tion and ii'-corrected model magnitude. The superscript 
0.1 means the rest- frame magnitude iiT-corrected to the 
redshift of 0.1, and will subsequently be dropped. We 
adopt the values of {g — i) color, concentration index 
(cin), and Petrosian radius i?pet computed for the galax- 
ies in KIAS-VAGC (Choi 2010, in prep.). The {g - i) 
color gradient was defined by the color difference be- 
tween the region with _Rgai < 0.5i?pct and the annulus 
with 0.5i?pct < ^gai < -Rpot, where i?gai is the galacto- 
centric radius and i?pot is the Petrosian radius estimated 
in the i-band image. To account for the effects of flatten- 
ing or inclination of galaxies, elliptical annuli were used 
to calculate the parameters. The (inverse) concentration 
index is defined by i?5o/i?9o, where R50 and Rgo are semi- 
major axis lengths of ellipses containing 50% and 90% of 
the Petrosian flux in the i-band image, respectively. 

The velocity dispersion values of galaxies are adopted 
from t he DR7 release of New York University (NYU) 
VAGG (jBlanton et al. 2005) and values of Ha equivalent 
width are taken from the DR7 release of Max-Planck- 
Institute for Astro physics (MPA)/John H opkins Univer- 
sity (JHU) VAGG (|Tremontiet al.l[200l . 

2.2. Satellite Systems of Isolated Galaxies 

Finding isolated satellite systems consists of two steps: 
identifying isolated galaxies and finding satellites asso- 
ciated with the isolated galaxies. We consider all SDSS 
galaxies with redshifts less than 0.03 as the host candi- 
dates. Only the relatively nearby galaxies are considered 
because the visual determination of the sense of rota- 
tion is feasible only when the galaxy image shows de- 
tailed internal features. For each galaxy its neighbors are 
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searched. The neighbors of a target galaxy with a r-band 
absolute magnitude Mr are the galaxies with an absolute 
magnitude brighter than M^+l and a velocity difference 
less than Av. We adopt Av = 1000 km s~^ for early-type 
target galaxies, and 750 km s~^ for late-type galaxies 
based on the distributions of pair wise velocity differe nce 
of the SP SS galaxies (see F ig. 2 of lPark fc Choill2009l and 
Fig. 1 of lWang et al.ll2009() . The neighbors are searched 
in the redshift interval from —0.005 to 0.034. When a 
target galaxy has a projected distance to its neighbor 
galaxies (i?noi) larger than max[rvir, target, ^^vir.nei] for aU 
neighbors, it is selected as an isolated gal axy. Here ryir 
is the virial radius of a galaxy (see eq. (5) of lPark fc Choil 
[20091) . 

Once isolated host candidates are determined, the 
satellites associated with them are searched for. A galaxy 
qualifies as a satellite if it has an absolute magnitude 
fainter than Mr{host)+l, velocity difference less than 
600 km s~^, and projected separation of the galaxy 



with the host R <min[i?nei ?-vir,host/(?'vir,host 



' vir,nciy 7 



i?,) 



i), 



i] where i?nei is a projected separation be- 



tween the host and host's neighbor. The first limit di- 
vides the distance between the host and neighbor in pro- 
portion to their virial radii. The second limit prohibits 
the galaxies within the virial radii of the neighbor galax- 
ies from being selected as satellites. In our SDSS galaxy 
sample we find 8904 satellites associated with 3515 iso- 
lated host galaxies. 

2.3. Spin Direction of Host Galaxies 

To determine the sense of rotation of host galaxies from 
their images we first assume that all spir al galaxies have 
trailing arms (jBinnev fc Tremainelll987l ) and determine 
the receding side of the major axis from their color im- 
ages. The part of minor axis behind the bulge often ap- 
pears smooth due to the foreground bulge stars while the 
side of the minor axis closer than the bulge often shows 
interstellar dust patches of relatively higher sharpness 
(jMihalas fc Binnevlll981f) . When it is not possible to de- 
termine the sense of rotation due to various reasons, we 
simply drop the host galaxy. Very early-type galaxies 
with smooth disks, very late-type galaxies with vanish- 
ing bulge, and nearly face-on or edge-on disk galaxies are 
difficult to assign the spin direction. After we remove all 
these cumbersome cases, we can determine the spin vec- 
tor directions for 165 host galaxies that have 535 satellite 
galaxies. 

We supplemented our sample by adding host galaxies 
of which rotation sense was determined fr om the spec- 
trosc opi c ob ser vation s in the literature (iRubin et al.l 
1981 [l985l [l99l iBro eils fc van Woerderi '1994' 
Han et al.lll995HTavlor et al . 1995; Rhcc & van Albada 
' 19971: iSo fuc 1997; Sofiie et al.l 
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Herbert-For t et al.l [2008; .Epinat et al.l l2008allbl : 
Swaters et al. 2009). We compiled 390 satellites asso- 
ciated with 116 hosts. There are 51 hosts in common 
between the hosts of which rotation sense was deter- 
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Fig. 1. — (a) Peculiar velocities of satellites as a function of 
absolute magnitudes of late-type host galaxies. Dashed-line en- 
velops indicate the velocity criteria for rejecting interlopers. (6) 
Distribution of absolute magnitudes of all hosts (open histogram) 
and of those having satellites within the velocity criteria (shaded 
histogram), (c) Distribution of peculiar velocities of all satellites 
(open histogram) and of those within the velocity criteria (shaded 
histogram). 

mined in this study and those in the literature. We 
have compared our determination of the rotation sense 
with that in the literature for these common objects, 
and found that there is only one galaxy (NGC 5894) 
that shows a disagreement between the two. We kept 
the rotation sense for this object determined from the 
spectroscopic observation in the literature. In total, 
we have 727 satellites assigned to 230 hosts whose spin 
vector directions are known. 

Among 727 satellites, we eliminated 45 spurious 
sources (e.g., faint fragments of bright galaxies, diffrac- 
tion spikes of bright stars) that were included in the list 
of satellites (originally included in the spectroscopic sam- 
ple of galaxies in SDSS) by inspecting color images of 
the satellites. In the result, 682 satellites with 221 host 
galaxies remain. 

It is noted that our list contains some satellite sys- 
tems whose host galaxies are early types because their 
spin directions were determined from the spectroscopic 
observations. Since the morphology of satellite galaxies 
is strongl y affected by the m orphology of host galaxies as 
found in lAnn et al.l ()2008[ ). we further divide the satel- 
lites into two categories (the cases of early- and late-type 
hosts) and obtain 197 late-type host galaxies having 579 
satellite galaxies and 24 early-type hosts with 103 satel- 
lites. Our analysis will be restricted to the satellite sys- 
tems with late-type hosts except Sill 

In Figure [1] we plot peculiar velocities of satellite 
galaxies as a function of absolute magnitudes of host 
galaxies. The peculiar velocity of a satellite is defined by, 
i^poc = 5'(wsat - I'host), where Wsat and Whost are observed 
velocities of the satellite and the host, respectively. S is 
defined by 



S 



-1, d < 90°, 
-1, otherwise. 
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Fig. 2. — Peculiar velocity distribution for all satellite galaxies 
associated with late- type host galaxies. The vertical dotted line 
indicates zero velocity. Dashed lines represent the double Gaussian 
fits, and the solid line indicates the sum of two Gaussian functions. 
The hatched histogram indicates the peculiar velocity distribution 
of satellite galaxies located near the major axis of the hosts {8 < 
45°). 

9 =1 PArcd — PAS | where PAred is a position an- 
gle of the receding side of the host's major axis, and 
PAS is the position angle of the satellite with respect to 
the host. All position angles are measured from north 
to east. Therefore, the satellites with i;pec > rotate 
around the host in the same direction as the rotation of 
the host (prograde orbit) and those with Wpcc < rotate 
in the opposite direction (retrograde orbit). 

Since the probability for a galaxy to become a satellite 
of a host decreases as the mass (luminosity) of the host 
decreases, we exclude some probable interlopers located 
outside the region bounded by dashed lines in Figure [TJ 
The dashed lines are determined as follows. The escape 
velocity of a satellite can be given by Wosc k (Mhost/-R)°'^ 
where Mhost is a mass of the host and i? is a projected 
distance to the host. The mass of the host galaxy is com- 
puted by Mhost = jL where 7 is the mass-to-light ratio 
and L is the r-band luminosity of the host. We assume 
that 7 (early) =27 (late) at the same r-band luminosity, 
and that 7 is consta nt with galaxy lum inosity for a given 
morphological type (|Park et al.l 120081) . Then we obtain 
VeBc{Mr) = uolO(-°-4/3)(*^'— *^o) j^ a giveu distance to 
the host galaxy (e.g., R — ^vir^host)- We adopt vq — 600 
km s~^ and Mq = — 21 mag, and plot Wosc versus Mr as 
dashed lines in Figure [TJ 

The virial radius of a host galaxy (rvir,host) is defined 
by the projected radius where the mean mass density p 
within the sphere of radius r is 200 times the critical 
density or 740 times the mean density of the universe, 
namely, 

rvir = (37i/47r/200pc)'/'. (2) 

Since we adopt n„, = 0.27, 200pc = 200p/rj,„ = 
740p where p is the mean density of the universe. 
The value of mean density of the universe, p = 
(0.0223 ± 0.0005)(7L)_2o(/i"^Mpc)-3, was adopted, 
where (7L)_2o is the mass of a late- type galaxy with 
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Fig. 3. — (a) Peculiar velocity ripoc vs. projected distance to 
the host galaxy. Median values in each radial bin and their errors 
are connected by solid lines. Peculiar velocity distribution for the 
satellites (b) at 0.1 < R/r^i,. host < 0-5 and (c) at 0.5 < R/r^i,. host- 
The vertical dotted lines indicate zero velocity. Solid lines indicate 
the best fit exponential functions of the distributions of satellites in 
prograde or retrograde orbits separately. Solid line for the satellites 
in prograde orbit at 0.1 < ^/''vir.host < 0-5 is the best fit Gaussian 
function. 



Mr = -20 (jPark et al.ll2008[ ). According to our formula 
the virial radii of galaxies with Mr = —19.5, —20.0, and 
—20.5 are 260, 300, and 350 h~^ kpc for early types, and 
210, 240, and 280 h^^ kpc for late types, respectively. 
By rejecting probable interlopers as shown in Figure [1] 
we finally secured 534 satellite galaxies associated with 
191 late-type hosts and 101 satellites associated with 24 
early-type hosts. 

We present the peculiar velocity distribution of all 534 
satellite galaxies in Figure [2] The fraction of galaxies in 
prograde orbit (/pro) is found to be 50 ±2%, and the me- 
dian value of the peculiar velocity distribution is — 1 ± 8 
km s~^. When we use only 262 satellites located near 
the major axis of the hosts {9 < 45°), /pro does not 
change and the median value of the distribution is ± 10 
km s^^. We fit the peculiar velocity distribution using 
the sum of two Gaussian functions, and find the best-fit 
function with mean values of (wpcc) — 7 km s~^(cr — 85 
km s"^) and (?;pec) = -22 km s-'^{a = 246 km s'^). 
When we fit the peculiar velocity distribution of satel- 
lites in prograde or retrograde orbits separately using a 
single Gaussian or exponential function, we find that the 
exponential function fits the distribution better than the 
Gaussian function. The reduced x^ for satellites in pro- 
grade orbit, are 2.1 and 0.5 for Gaussian and exponential 
functions, respectively. The reduced x^ for satellites in 
retrograde orbit are 2.5 and 0.7 for Gaussian and expo- 
nential functions, respectively. 

To study whether or not the peculiar velocity distribu- 
tion changes depending on the distance to the host we 
show the distribution for satellites at 0.1 < R/r^ir host < 
0.5 or at 0.5 < i?/rvir,host separately in Figure [3| We 
use Gaussian and exponential functions to fit the distri- 
butions of satellites in prograde and retrograde orbits, 
and show best-fit functions in Figure [S^b-c). It is seen 
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Fig. 4. — Absolute magnitude vs. redshift for late- type host 
galaxies. The dashed line is the apparent magnitude limit {rrir = 
17.77) computed from eq. l(lj in Choi ct al. (20071). The satellites 
associated with the host galaxies within the solid box are selected 
as a luminosity subsample of satellites. 
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Fig. 5. — Absolute magnitude vs. redshift for early- type satel- 
lites in prograde orbit (Ep), for early- type satellites in retrograde 
orbit (Er), for late- type satellites in prograde orbit (Lp), and for 
late- type satellites in retrograde orbit (Lr). Satellites of the host 
galaxies in the luminosity subsample (—18.5 > M^ ij^gt > —21.5) 
are filled circles, while the rests are crosses. Number outside the 
parenthesis means the number of galaxies in each panel, and the 
number in the parenthesis indicates the number of galaxies within 
the luminosity subsample. 

that the best-fit function for satellites in retrograde or- 
bit is exponential both in the inner and outer regions, 
while that for satellites in prograde orbit changes from 
Gaussian (inner region) to exponential (outer region). 

We plot r-band absolute magnitudes of host galaxies 
and their satellites as a function of redshift in Figures 
13] and We made a luminosity subsample of satellites 
whose hosts have absolute magnitudes in the range of 
-18.5 > Mr > -21.5. 

In Figure [U we plot peculiar velocities of satellites as 
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Fig. 6. — Peculiar velocity iipcc vs. projected distance to the host 
galaxy for early- type satellites in prograde orbit (Ep), for early- 
type satellites in retrograde orbit (Er), for late- type satellites in 
prograde orbit (Lp), and for late-type satellites in retrograde orbit 
(Lr) in the luminosity subsample. Global median values are shown 
by horizontal dot-dashed lines, and the median value in each radial 
bin is connected by the solid line. Vertical dashed line indicates 
one virial radius of the host. 

a function of a distance to host galaxies. It shows that 
global median values of iJpoc do not significantly vary 
depending on the orbit or morphology. However, there 
exists an interesting trend that f poc of the late- type satel- 
lites in prograde orbit (Lp) increases as they approach 
the host, while Wpoc of those in retrograde orbit (Lr) de- 
creases. At the merger scales with i?/rvir,host < 0.1 the 
orbital velocity decreases for both Lp and Lr cases. 

3. RESULTS 
3.1. Spatial Distribution of Satellites 

In Figure [71 we plot the radial distribution of satellites 
in prograde and retrograde orbits. We also show the 
fraction of galaxies in prograde orbit (/pro) as a function 
of the projected distance to host galaxies. The top panel 
shows that satellites in prograde orbits are located closer 
to their hosts compared to those in retrograde orbits for 
the case of early-type satellites. However, we do not find 
such dependence of the spatial distribution of late-type 
satellites on their orbits (see panel c). As a result, /pro 
of early types increases as R decreases, but /pro of late 
types remains constant (~ 0.5) within the error bar. 

To study whether or not satellites are concentrated on 
the plane of host galaxies' disk, we plot the ratio of the 
dispersion of satellite distribution along {a{cos9)) and 
perpendicular (a{sm6)) to the major axis of the host as 
a function of the ratio of minor axis to the major axis of 
the host {{b/a)iioST) in Figured It shows neither of the 
distributions of prograde or retrograde orbit satellites is 
aligned with late-type host galaxies. 

3.2. Luminosity 

Figure [HI shows absolute magnitudes of the satellites in 
the luminosity subsample as a function of the projected 
distance to host galaxies. Each panel represents differ- 
ent subsamples depending on morphology and orbit. For 
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Fig. 7. — Distribution of early-type satellites in prograde (filled 
triangles) and retrograde (open upside-down triangles) orbits (a) 
and the number fraction of satellites in prograde orbit (/pro) among 
early-type satellites (b) as a function of the projected distance to 
the host galaxies. Distribution of late-type satellites in prograde 
(filled triangles) and retrograde (open upside-down triangles) orbits 
(c) and /pro among late- type satellites (d). 
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Fig. 8. — Ratio of the dispersions of satellite distribution along 
((t(cos9)) and perpendicular {a-{sin9)) to the major axis of hosts as 
a function of the host axis ratio. Solid and dashed lines indicate 
satellites in prograde and retrograde orbits, respectively. 



both early and late types, galaxies in retrograde orbit, 
on average, appear to be slightly brighter than those in 
prograde orbit. A statistical test with Monte Carlo simu- 
lation to compute the significance of the difference in the 
median of the physical quantity (to be discussed in the 
end of this section) confirms this difference for only late 
types with a significance level of 97%. When we compare 
satellites in each radial bin, the Monte Carlo test indi- 
cates with a significance level of 99% that median abso- 
lute magnitudes of late-type galaxies in the outer region 
(R > Tvir.host) are different depending on the orbit. 

To investigate the difference in the galaxy luminosity 
depending on the orbit in more detail, we show the dis- 
tribution of absolute magnitudes in Figure [TOl We count 
the number of satellites within the redshift limit defined 
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Fig. 9. — Same as Fig. [G] but for absolute magnitude Mr 
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Fig. 10. — Distribution of absolute magnitudes for all satellite 
galaxies. Solid and dotted lines, respectively, indicate luminosity 
functions (LFs) for late-type satellites in prograde and retrograde 
orbits. Dashed and dot-dashed lines denote LFs for early-type 
satellites in prograde and retrograde orbits, respectively. 

by the apparent magnitude limit in Figure [5l and do not 
correct the count for different volumes. It is seen that 
bright satellites (—18 > Mr > —20) in retrograde orbit 
are more numerous compared to those in prograde or- 
bit. The test with Monte Carlo simulation confirms this 
result with a significance level of 97% in the sense that 
median absolute magnitudes of late types in prograde 
and retrograde orbits are significantly different. 

We summarize the results of statistical comparison of 
several physical parameters between satellites in pro- 
grade and retrograde orbits and between satellites in the 
inner and outer regions in Tables [1] and [21 respectively. 

The statistical significance of the difference in the me- 
dian or the dispersion of physical quantities, is computed 
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Fig. 11. — Same as Fig. |6l but for (u — r) color. 

by Monte Carlo simulations. For a given physical pa- 
rameter, we construct two subsamples (with the same 
number of galaxies as in the actual prograde and retro- 
grade samples) by randomly drawing the observed pa- 
rameters from the whole early-type or late-type satellite 
galaxy sample and compute the median and the disper- 
sion of each subsamplc. The resulting two subsamples 
will have the same medians and standard deviations on 
average. These random subsamples will tell us about 
whether or not the difference in the parameter between 
prograde and retrograde satellite samples is statistically 
significant. We construct 1000 trial data sets and com- 
pute the fraction of simulated data sets in which the dif- 
ference of the physical quantity is larger than or equal to 
that based on the real data (fsim>obs)- The significance 
levels of the difference defined by 100(1— fsim>obs) (%) 
are given in Tables. We have applied the statistical tests 
to all the physical parameters, but only show statisti- 
cally significant results (significance level above 95%) in 
the Tables. 

3.3. Star Formation Activity Parameters 

Figure [TT] shows {u — r) colors of satellite galaxies di- 
vided by their morphologies and orbital motion as a func- 
tion of the projected distance to the hosts. The late-type 
satellites in prograde and retrograde orbits have a similar 
global median value, and both show a hint of an increase 
in {u — r) color with decreasing R. In the Monte Carlo 
test for the radial gradient we divide each sample into 
inner (i? < Tvir.host) and outer {R > rvir,host) regions, 
and test if the difference in the physical parameter in 
two regions is statistically significant. The random sub- 
samples are generated in such a way that a parameter 
value is randomly drawn from the combined sample and 
assigned to a galaxy in inner or outer region. The re- 
sulting subsamples will have no radial gradient in the 
parameter and have the same dispersions on avergae. A 
test with such Monte Carlo simulation confirms that the 
change of (u — r) color with R for late types in prograde 
orbit (see Tab. ^ is significant at a 99% confidence level. 
It is noted that the increase of {u — r) colors for those 
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Fig. 12. — Same as Fig. [G] but for equivalent width of Hq line. 
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Fig. 13. — Same as Fig. (6] but for (g — i) color gradient. 

in retrograde orbit is very strong in the innermost region 
(i? < O.lTvir.host) even though the overall gradient is not 
signifiant. 

In Figure I12[ we plot the equivalent width of Ha lines 
of satellite galaxies, which is a useful measure of star for- 
mation activity (SFA). For late- type satellites, the me- 
dian values of W(Ha) for two subsamples are similar, and 
those do not seem to change with R for both samples even 
though their {u — r) colors become redder at smaller R as 
seen in Figure [TT] However, the test with Monte Carlo 
simulation shows that the difference in the dispersion of 
W(Ha) between late-type satellites in prograde and ret- 
rograde orbits is remarkable with a significance level of 
of 97% (see Tab. [T]). The difference in the dispersion 
of W(Ha) of retrograde orbit satellites in the inner and 
outer regions is also noted (see Tab. [2]). 

Figure [13] represents [g — i) color gradients of satel- 
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TABLE 1 

Statistical Comparison between Satellites in Prograde and Retrograde Orbits 



Parameter Type MC,, 



MC, 



d.all 



(1) 



(2) 



(3) 



(4) 



MC„,, 
(5) 



MCd,i 
(6) 



MCn 



MCj Qut Remarks 



(7) 



(8) 



(9) 



Mr 


L 


97 


58 


49 


40 


W(Ha) 


L 


56 


49 


33 


11 


A(3 - i) 


L 


95 


98 


66 


72 


Cin 


L 


38 


98 


47 


88 



99 



38 



69 


Fig. [9] 


97 


Fig.[T2] 


98 


Fig.[T3] 


97 


Fig.[Tl 



Column descriptions. (1): Physical parameters. (2): Morphological types (E: early types, L: late types). (3); Significance level (%) of the 
difference of the median of the physical parameter between satellites in prograde and retrograde orbits determined from Monte Carlo 
simulation. (4): Significance level (%) of the difference of the dispersion of the physical parameter between satellites in prograde and 
retrograde orbits determined from Monte Carlo simulation. (5): Same as column (3), but for satellites in the inner region {R < Irvir host)- 
(6); Same as column (4), but for satellites in the inner region. (7): Same as column (3), but for satellites in the inner region. (8): Same as 
column (4), but for satellites in the outer region (/? > Irvir host)- (9)- Relevant figure. Statistically significant values (> 95%) are 
represented in bold face. 



TABLE 2 
Statistical Comparison between Satellites in the inner 

AND outer regions 



Parameter 


Type 


-LVlv_'ni,pro 


MCd,pro 


MCm,ret 


MCd,rct 


Remarks? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


{u-r) 


L 


99 


67 


88 


65 


Fig-in 


W(Ha) 


L 


13 


12 


64 


97 


Fig-E] 


A(5 - i) 
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96 


61 


4 


43 


Fig. [13] 


Rpct 
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97 


77 
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24 


Fig. [H 


{u — r) 
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99 


98 


11 


27 


Fig. El 


{u — r) 
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93 


98 
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67 


Fig. [18] 



Column descriptions. (1): Physical parameters. (2): 
Morphological types (E: early types, L: late types). (3); 
Significance level (%) of the difference of the median of the 
physical parameter of prograde orbit satellites in the inner and 
outer regions determined from Monte Carlo simulation. (4): 
Significance level (%) of the difference of the dispersion of the 
physical parameter of prograde orbit satellites in the inner and 
outer regions determined from Monte Carlo simulation. (5): 
Same as column (3), but for retrograde orbit satellites. (6): Same 
as column (4), but for retrograde orbit satellites. (7): Relevant 
figure. Statistically significant values (> 95%) are represented in 
bold face. 

lite galaxies. The positive {g — i) color gradient means a 
bluer color in the central region of a galaxy than that in 
the outer region. For early types, the difference of me- 
dian {g — i) color gradients of satellites in prograde orbit 
(Ep) between in the inner and outer regions is notice- 
able, which is supported by the test with Monte Carlo 
simulation (see Tab. [1]). For late types, the test with 
Monte Carlo simulation suggests that the median and 
its dispersion of {g — i) color gradients of satellites in 
prograde and retrograde orbits are different significantly 
in the whole/outer region. 

3.4. Structure Parameters 

Figure [14] shows the distribution of the concentration 
index Cjn for satellite galaxies divided by their morpholo- 
gies and orbital motion. Small value of Cjn means that 
the light of a galaxy is more concentrated in the cen- 
tral region. It is found that the global median and the 
scatter of cin values for early-type satellites are not sig- 
nificantly different depending on their orbits. However, 
it is noted that there are few early types in retrograde 
orbit with large Cjn values (> 0.4) in the inner region 
{R < 0.2rvir,host), and few early types in prograde orbit 
with large c-m values (> 0.4) are found in the outer region 
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Fig. 14. — Same as Fig. [6] but for concentration index ci^. 

{R > r-vir,host)- For late types, it is seen that the differ- 
ence in the dispersion of c-m between late-type satellites 
in the prograde and retrograde orbits (see Tab. [ij. How- 
ever, the global median values of Cm for two subsamples 
are similar, and the median values for both subsamples 
do not show any significant variation with R. 

In Figure (TSl we plot central velocity dispersions of 
satellite galaxies. The median values for early types in 
prograde and retrograde orbits are similar and do not 
show any significant variation with R. For late types, 
any significant difference between two subsamples of late- 
type satellites is not found except for the scatter of a 
values. 

Figure \W\ represents the i-band Petrosian radius i?pot 
of satellite galaxies. The i?pot for the early-type satellites 
in the inner region {R < r-vir.host) appears to be larger 
than that in the outer region (i? > rvir,host) (see Tab. 
[2|), while no significant difference between the late- type 
satellites in prograde and retrograde orbits is not seen. 



4. DISCUSSION 



4.1. 



Different Galaxy Samples 

We found in Figures fl2l[T3] and Table [1] that the median 
or its dispersion of W(Ha) and {g — i) color gradients can 
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Fig. 15. — Same as Fig. |6l but for central velocity dispersion a. 
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Fig. 16. — Same as Fig. |6l but for Petrosian radius ijpot- 

be different depending on the orbital motion of late-type 
satellites. The SFA of early- type satellites in retrograde 
orbit (Er) also appears to be stronger than those in pro- 
grade orbit in the outer region (rvir,host < ^) as seen in 
Figures [Til and fT2l Moreover, (u — r) colors of late-type 
satellites increase as they approach their host galaxies 
for both satellites in prograde and retrograde orbits, and 
(u — r) colors of late-type satellites in retrograde orbit 
in the innermost region (J? < O.lrvir.host) are remarkably 
large. These differences depending on the orbit can be 
attributed to the different origin of satellite galaxies, or 
to the different strength of physical processes that they 
experience through interactions with host galaxies. 

To check whether or not these results are strongly af- 
fected by the satellites associated to the host galaxies 
whose spin direction is determined from the SDSS color 
images in this study, in Figure [171 we plot (u — r) col- 



y.u 


: Ep 


' _' ': 


: Er 







2.5 


- 


! - 


- 


• 


: 


2.0 


■ 




- 




- 


*. 












1.5 


- 




- 






1.0 








1 


: 


y.o 


: Lp 




: Lr 






2.5 


. 


\ - 


- t 




2.0 


- 


\ _ - 


- 


\" •• :i; 




1 5 


■ . .i-=-Tr>^ 


~S4^'- ^ - 


.._._. 


\t* •''■^ib* 


-^?" 






^c-- : 




1.0 


, 




- 




" 



0.1 



1.0 



0.1 



1.0 



R/fvi 



R/r„, 



Fig. 17. — Same as Fig. 1111 but for the satellites associated to 
host galaxies whose spin direction is determined from the spectro- 
scopic observations. 
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Fig. 18. — Same as Fig. 1111 but for the satellites associated to 
early-type hosts. 

ors versus the projected distance to the hosts using the 
satellites associated to the host galaxies whose spin direc- 
tion is determined from the spectroscopic observations. 
For both early and late types, we do not find any sig- 
nificant difference of (u — r) colors between satellites in 
prograde and retrograde orbits, which is consistent with 
the result when we use the combined sample of satellites 
in the host galaxies whose spin direction is determined 
from the color images or the spectroscopic observations 
(see Fig. [TT|) . In addition, the late- type satellites in pro- 
grade orbits (Lp) show an increase in {u — r) color as 
decreasing R, which is again similar to the case of late 
types in prograde orbits using the combined sample of 
satellites (see Lp in Fig. [TTJand Tab. [5]). In the result, 
the results of SFA using the combined sample of satel- 
lites are not affected strongly by the satellites associated 
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to the host galaxies whose spin direction is determined 
from the color images in this study. 

Since the physical parameters of satellite galaxies can 
be stron gly affected by t he morphology of host galax- 
ies (e.g.. lAnn e t al." 2008"). it is also interesting to com- 
pare satellite systems that have early- and late-type host 
galaxies. To investigate the SFA of satellites around 
early-type hosts we use a sample of 101 satellites asso- 
ciated with 24 early-type hosts, and show their {u — r) 
colors in Figure UHl It is seen that (m — r) colors of late- 
type satellites become redder as R decreases, which was 
also seen in Figure [Til or [TTl For early-type satellites, the 
small number of satellites prevents us from drawing any 
firm conclusion. 

4.2. Comparison with Other Galaxy Systems 

Previously, iParket al.' (2008) and lPark fc Choil ()2009l ) 
found that galaxy properties strongly depend on the dis- 
tance and the morphology of the nearest neighbor galaxy. 
For example, the SFA of galaxies is enhanced as they 
approach late-type neighbors of comparable luminosity, 
while it is reduced as they approach early-type neighbors. 
The structure parameters of galaxies change with the 
distance to the nearest neighbor galaxy in the same di- 
rection independently of neighbor's morphology. Galaxy 
properties start to change when galaxies approach their 
nearest neighbor galaxy closer than the virial radius 
(rvir,nci) of the nearest neighbor galaxy and to change 
significantly at i? < 0.05rvir,nci where the galaxies in pair 
start to merge. They interpreted this phenomenon as a 
result of hydrodynamic interactions within the virial ra- 
dius of the galaxy plus dark halo system. 

The dependence of satellite properties as a function 
of the distance to the host galaxy in this study can 
be compared with the case of late-type neighbors in 
iPark fc Choil ()2009[ ). Since we restrict our analysis to 
the case of late-type host galaxy, it is expected that the 
SFA of both early- and late-type satellites increases as 
they approach the host. However, we found in i )3.3l that 
(u—r) color becomes redder (Fig. [TT]) and W{Ha) hardly 
changes (Fig. [T2|) even when late-type satellites approach 
their late-type hosts. This may be becau se, unlike the 
galaxies considered bv lPark fc Choil ()2009f ). the satellites 
in this study are much fainter and less massive than their 
hosts. As a result they may not be able to pick up cold 
gas from their host much, and the enhancement of SFA 
through interactions with their host is counterbalanced 
by the quenching by the hot halo gas of the host. 

On the other hand, the environment of satellite galax- 
ies can be compared with that of galaxies in galaxy clus- 
ters because the member galaxies in a galaxy cluster are 
bounded within the gravitational potential well of the 
cluster as the s atellites are bounded by their host galaxy. 
iPark &: Hwana (2009. ) studied the dependence of galaxy 
properties on the clustercentric radius and the environ- 
ment attributed to the nearest neighbor galaxy using the 
SDSS galaxies associated with the Abell galaxy clusters. 
They found that the SFA of late-type galaxies decreases 
as they approach the cluster center, and suggested that 
hydrodynamic interactions with nearby early-type galax- 
ies in clusters is the main drive to quenching the SFA of 
late types. Most cluster galaxies are early types includ- 
ing the brightest cluster galaxies (or cD galaxies), which 
might be comparable with the case of early-type host 



galaxy in the satellite system as seen in Figure [TH] The 
late-type satellites (bottom panels in Fig. \TE\\ become 
redder as they approach early-type host galaxies, which 
agrees to the case of gala xy clusters. 

Previously, iZinnI (|1993l ) pointed out that there are two 
distinct populations in the halo globular clusters in the 
Milky Way: old and young halo clusters. Two popula- 
tions were originally divided by the horizontal branch 
morphology, but were found to be also separated by 
their orbital motion in the Milky Way: progradc or- 
bit for old halo clusters and retrograde orbit for young 
halo clusters (I Zinni i 19931: Ivan den Berghlll993l : iLee et al.l 
l2007f) . Therefore, the formation history of old halo clus- 
ters with prograde orbit is thought to be different from 
that of young hal o cl usters with retrograd e orbit (e.g., 
lEggen et all 119621 vs. iSearle fc ZImI 119781 ). For exam- 
ple, it was suggested that the clusters in retrograde or- 
bit may have bee n accreted into the Milky Way (e.g., 
iDinescu et al.llT999 !). and further that these clusters and 
Galactic dSphs may have the same origin in accretion or 
in partially disru ption of formerly large paren t satellite 
galaxies (e.g.. Ivan den Bergh fc M ackev '2004). There- 
fore, two populations of satellite galaxies with different 
orbits in this study are also expected to have formed and 
evolved in different paths, which may be responsible for 
the different SFA depending on the orbit. The Galactic 
satellite galaxies also can give us important clues for test- 
ing this scenario, but it is not achievable at this moment 
because of small number (~ 10) of Galactic satellites 
whose orbital motion is determined. 

5. CONCLUSIONS 

We have studied the orbital dependence of various 
galaxy properties in the satellite systems of galaxies in 
SDSS. Our primary results are summarized as follows. 

1 . We have found the satellite systems of galaxies at 
z < 0.03, which consist of 8904 satellites associated 
with 3515 isolated host galaxies. 

2. We have determined the spin direction of late- type 
host galaxies by inspecting the SDSS color images. 
By combining our results with those from the pre- 
vious spectroscopic observations in the literature, 
we obtained 635 satellites associated with 215 host 
galaxies whose spin directions are known. 

3. The number ratio of satellites in prograde orbit to 
those in retrograde orbit is nearly equal to one, 
and the peculiar velocity distribution is found to 
be symmetric. The number ratio of satellites in 
prograde orbit to those in retrograde orbit (/pro) 
appears to increase as the projected distance to 
the host galaxy decreases. However, we do not find 
any significant difference in the azimuthal distri- 
bution of satellites between those in prograde and 
retrograde orbits. 

4. The satellites in retrograde orbit, on average, ap- 
pear to be slightly brighter than those in prograde 
orbit. The latter seems to get brighter as they ap- 
proach late-type hosts, while the former does not 
show any significant change with the distance to 
the hosts. 
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5. For early- type satellites in retrograde orbit, the 
SFA seems to decrease as they approach host galax- 
ies. For late-type satellites, {u — r) colors of both 
samples in prograde and retrograde orbits appear 
to increase as they approach the hosts, while the 
radial variations of the equivalent width of Ha lines 
for both samples are not significant. 

6. The structure parameters of satellites do not seem 
to be different significantly between galaxies in pro- 
grade and retrograde orbits. 
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